arXiv: 1502.00586v2 [hep-ex] 20 May 2015 


EUROPEAN ORGANISATION FOR NUCLEAR RESEARCH (CERN) 



A 

T 

L 

A 

S 



CERN-PH-EP-2014-284 

Submitted to: Physical Review D 


Observation of top-quark pair production in association with a 
photon and measurement of the ££7 production cross section in 
pp collisions at Vs = 7 TeV using the ATLAS detector 


The ATLAS Collaboration 


Abstract 


A search is performed for top-quark pairs (tt) produced together with a photon (7) with trans¬ 
verse energy greater than 20 GeV using a sample of tt candidate events in final states with jets, 
missing transverse momentum, and one isolated electron or muon. The dataset used corresponds 
to an integrated luminosity of 4.59 fb _1 of proton-proton collisions at a center-of-mass energy of 
7 TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. In total 140 and 222 
tt') candidate events are observed in the electron and muon channels, to be compared to the ex¬ 
pectation of 79 ± 26 and 120 ± 39 non -tlj background events respectively. The production of 
tt') events is observed with a significance of 5.3 standard deviations away from the null hypoth¬ 
esis. The tt') production cross section times the branching ratio (BR) of the single-lepton decay 
channel is measured in a fiducial kinematic region within the ATLAS acceptance. The measured 
value is < 7 ^ x BR = 63 ± 8 (stat.) ± J 3 (syst.) ± 1 (lumi.) ft per lepton flavor, in good agreement with 
the leading-order theoretical calculation normalized to the next-to-leading-order theoretical prediction 
Of 48 ± 10 ft. 
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Observation of top-quark pair production in association with a photon and 
measurement of the tt'y production cross section in pp collisions at y^s = 7 TeV using 

the ATLAS detector 

The ATLAS Collaboration 

A search is performed for top-quark pairs ( tt) produced together with a photon ( 7 ) with transverse 
energy greater than 20 GeV using a sample of tt candidate events in final states with jets, missing 
transverse momentum, and one isolated electron or muon. The dataset used corresponds to an 
integrated luminosity of 4.59 fb -1 of proton-proton collisions at a center-of-mass energy of 7 TeV 
recorded by the ATLAS detector at the CERN Large Hadron Collider. In total 140 and 222 tt'y 
candidate events are observed in the electron and muon channels, to be compared to the expectation 
of 79 ± 26 and 120 ± 39 non-tty background events respectively. The production of tty events 
is observed with a significance of 5.3 standard deviations away from the null hypothesis. The 
tty production cross section times the branching ratio (BR) of the single-lepton decay channel is 
measured in a fiducial kinematic region within the ATLAS acceptance. The measured value is 
afj*. x BR = 63 ± 8 (stat.) jjjg (syst.) ± 1 (lumi.) fb per lepton flavor, in good agreement with the 
leading-order theoretical calculation normalized to the next-to-leading-order theoretical prediction 
of 48 ± 10 fb. 

PACS numbers: 14.65.Ha, 12.60.Jv, 13.85.Qk, 14.80.Ly 


I. INTRODUCTION 

Due to its large mass, the top-quark is speculated to 
play a special role in electroweak symmetry breaking 
(EWSB). New physics connected with EWSB can man¬ 
ifest itself in top-quark observables. For instance, top- 
quark couplings can be modified significantly in some ex¬ 
tensions of the Standard Model (SM). A measured yield 
of top-quark pair production in association with a pho¬ 
ton (tty) can constrain models of new physics, for exam¬ 
ple those with composite top-quarks [1], or with excited 
top-quark production, followed by the radiative decay 
t* ty. The tty coupling may be determined via an 
analysis of direct production of top-quark pairs in asso¬ 
ciation with a photon, evidence of which was first re¬ 
ported [2] by the CDF collaboration. 

In this paper, observation of top-quark pair produc¬ 
tion in association with a photon in proton-proton (pp) 
collisions at a center-of-mass energy of s/s = 7 TeV is 
presented using the full 2011 ATLAS data sample, which 
corresponds to an integrated luminosity of 4.59 fb -1 . 
This analysis is performed on tt candidate events in the 
lepton plus jets final state. The tty candidates are the 
subset of tt candidate events with an additional pho¬ 
ton. The measurement of the tty production cross section 
times the branching ratio (BR) of the single-lepton decay 
channel (iv^qq’bby, where £ is an electron or muon) is re¬ 
ported in a fiducial kinematic region within the ATLAS 
acceptance. 

The paper is organized as follows. The ATLAS detec¬ 
tor is briefly described in Sec. II. The data and Monte 
Carlo simulation samples used in the analysis are de¬ 
scribed in Sec. Ill, followed by a description of the event 
selection in Sec. IV. The definition of the fiducial phase 
space used in the measurement is presented in Sec. V. 
The cross section is extracted from a template-based pro¬ 


file likelihood fit using the photon track-isolation distri¬ 
bution as the discriminating variable. Section VI details 
the overall strategy of the measurement, and describes 
how prompt-photon and background templates are ob¬ 
tained. Background estimates are discussed in Sec. VII. 
An overview of the systematic uncertainties in the mea¬ 
surement is presented in Sec. VIII. Section IX presents 
the results of the measurement, followed by conclusions 
in Sec. X. 


II. DETECTOR 

A detailed description of the ATLAS detector can be 
found in Ref. [3] . The innermost part of the detector is a 
tracking system that is immersed in a 2 T axial magnetic 
field and measures the momentum of charged particles 
within a pseudorapidity range of |?y| < 2.5 a . The inner 
detector (ID) comprises silicon pixel and microstrip de¬ 
tectors, and a transition radiation tracker. The calorime¬ 
ter system is composed of sampling electromagnetic and 
hadronic compartments with either liquid argon or scin¬ 
tillator tiles as the active media. It resides outside the 
ID, covering \r/\ < 4.9. The outermost system is a muon 
spectrometer that is used to identify and measure the 
momentum of muons in a toroidal magnetic field in the 
region \rj\ < 2.7, with detectors used for triggering within 
\r]\ < 2.4. A three-level trigger system selects the po¬ 
tentially interesting events that are recorded for offline 
analysis. 


a ATLAS uses a right-handed coordinate system with its origin 
at the nominal interaction point (IP) in the center of the de¬ 
tector and the z-axis along the beam pipe. The x-axis points 
from the IP to the center of the LHC ring, and the y-axis points 
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III. DATA AND MONTE CARLO SAMPLES 

Data recorded by the ATLAS detector in 2011 in pp 
collisions at y/s = 7 TeV are considered for analysis. Re¬ 
quirements are imposed on the collected data to ensure 
the quality of the beam conditions and detector perfor¬ 
mance. The total integrated luminosity of the analyzed 
data sample is C = 4.59 ± 0.08 fb -1 [4]. 

Monte Carlo simulation samples are used to study sig¬ 
nal and background processes, using the ATLAS detector 
simulation [5] based on the GEANT4 program [ 6 ] . To simu¬ 
late effects of multiple pp interactions per bunch crossing 
(‘pile-up’), all Monte Carlo events are overlaid with ad¬ 
ditional inelastic events generated with PYTHIA [7] using 
the AMBT1 set of parameters (tune) [ 8 ]. The events are 
then reweighted to match the distribution of the mean 
number of interactions per bunch crossing in the data. 
Simulated events are reconstructed in the same manner 
as the data. 

Signal tty events with single-lepton {tv^qq'bb^, 
£ = e, /x, t) or dilepton (£u(£'upbir/, £/£' = e, /z, r) final 
states are simulated with two independent leading-order 
(LO) matrix element (ME) Monte Carlo generators, 
WHIZARD vl.93 [9, 10] and MadGraph v5.1.5.12 [11], 
both using the CTEQ6L1 [12] LO parton distribution 
function (PDF) set. Both calculations take into account 
interference effects between radiative top-quark produc¬ 
tion and decay processes. Details on the generator-level 
settings of the two signal Monte Carlo samples are 
available in Sec. Al. In the tty and inclusive ti samples 
the top-quark mass is set to m t = 172.5 GeV. 

The WHIZARD sample is interfaced to 
HERWIG v6.520 [13] for the parton showering and 
JIMMY 4.31 [14] is used for the underlying-event sim¬ 
ulation. The AUET2 tune [15] is used. The MadGraph 
sample is interfaced to either the PYTHIA v6.425 
parton shower using the PERUGIA 2011 C tune [16], or 
with HERWIG v6.520 and JIMMY 4.31 for the parton 
showering and the underlying-event simulations respec¬ 
tively. PYTHIA QED final-state radiation (FSR) from 
charged hadrons and leptons is switched off and instead 
PHOTOS v2.15 [17] is used. 

To compare with the experimental measurement, the 
LO calculations of WHIZARD and MadGraph are normal¬ 
ized to the next-to-leading-order (NLO) cross section, 
obtained for y/s = 7 TeV at the renormalization and 
factorization scales of mt . The NLO QCD calculation of 
top-quark pair production in association with a hard pho¬ 
ton is detailed in Sec. A 2. The systematic uncertainty on 
the NLO cross section is obtained by simultaneous renor¬ 
malization and factorization scale variations by a factor 


upward. Cylindrical coordinates (r, 0) are used in the trans¬ 
verse plane, 0 being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 6 as 
i] = — lntan(0/2). Transverse momentum and energy are de¬ 
fined as pry = p sin 0 and E-y = E sin 0. respectively. 


of two (mt /2 and 2 mt) around the central value ( mt ), 
and is calculated to be 20% [18]. The NLO/LO correc¬ 
tion (AT-factor) calculation is performed in a phase-space 
region close to the one defined by the analysis kinematic 
selection criteria (see Sec. A2 for details). The depen¬ 
dence of the A'-factor on the kinematic variables is small 
compared to the scale uncertainty [18]. 

The effect of the variations of photon radiation settings 
in MadGraph is studied using a sample generated with a 
minimum separation in p-cft space between the photon 
and any other particle of A R > 0.05 b instead of A R > 
0.2 used in the default sample (see Sec. Al). For this 
sample, PYTHIA QED FSR is switched off and no addi¬ 
tional photon radiation is produced by PHOTOS v2.15. In 
addition to the default MadGraph+PYTHIA Monte Carlo 
sample generated at the scale of m*, samples at scales of 
m t /2 and 2 m t are produced to study the effect of scale 
variations. 

The simulated sample for inclusive ti production 
is generated with MC0NL0 v3.1 [19, 20] (NLO ME 
2 —» 2) interfaced to HERWIG v6.520 for the parton 
showering and fragmentation and to JIMMY 4.31 for 
underlying-event simulation, using the CTEQ 6.6 [21] 
PDF set, with additional photon radiation simulated 
with PHOTOS v2.15. This sample is used to validate dis¬ 
tributions of kinematic variables in ti candidate events 
as described in Sec. IV. 

Initial- and final-state QCD radiation (ISR/FSR) vari¬ 
ations are studied using inclusive ti samples generated 
with AcerMC v3.8 [22] interfaced to PYTHIA v6.425 with 
the CTEQ6L1 PDF set. In these samples the pa¬ 
rameters that control the amount of ISR/FSR are set 
to values consistent with the PERUGIA Hard/Soft tune 
in a range given by current experimental data [23]. 
AcerMC v3.8 ti samples interfaced to PYTHIA v6.425 
are also used to study variations of color reconnection 
using the PERUGIA 2011 C and PERUGIA 2011 NO CR 
tunes [16]. The underlying-event variations are stud¬ 
ied using AcerMC v3.8 interfaced to PYTHIA v6.425 
with two different underlying-event settings of the 
AUET2B [24] PYTHIA generator tune. In all these 
AcerMC v3.8 samples, photon radiation is simulated 
with PHOTOS v2.15 [17]. The inclusive ti signal sam¬ 
ples are normalized to a predicted Standard Model ti 
cross section of a t t = 177l 1 b pb for a top-quark mass of 
172.5 GeV, as obtained at next-to-next-to-leading or¬ 
der (NNLO) in QCD including resummation of next- 
to-next-to-leading-logarithmic (NNLL) soft gluon terms 
with Top++ v2.0 [25-30]. 

Background samples of W and Z bosons (includ¬ 
ing W + bb and Z + bb processes) are generated with 
ALPGEN v2.13 [31] interfaced to HERWIG v6.520, using 
the CTEQ6L1 PDF set. The ALPGEN matrix elements 
include diagrams with up to five additional partons. The 


b A R = <J (Ay ) 2 -p (A?;) 2 , where A 77 (Aq) is the separation in // 
(0) between the objects in the //- 0 space. 
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MLM [31] parton-jet matching scheme is applied to avoid 
double counting of configurations generated by both the 
parton shower and the LO matrix-element calculation. 
In addition, overlap between heavy-flavor quarks that 
originate from ME production and those that originate 
from the parton shower is removed. Diboson (WW, WZ , 
and ZZ) production is modeled using HERWIG v6.520 
and the MRST LO** PDF set [32], The IVy+jets and 
A 7 +jets (with up to three partons including 66 , cc, c) 
processes are generated with SHERPA v 1.4.0 [33] and 
the CT10 [34] NLO PDF set. Single-top-quark pro¬ 
duction is modeled using AcerMC in the f-channel and 
MC0NL0 v3.41 [35] for the Wt- and [36] s-cliannels. 

Multijet samples with jet px thresholds of 17, 35 and 
70 GeV are generated using PYTHIA v6.421 with the 
AUET2B [24] generator tune. 


IV. OBJECT AND EVENT SELECTION 

Events for the analysis are selected by requiring a high- 
px single-electron or single-muon trigger [37] for the elec¬ 
tron and muon channels respectively. The px threshold 
for the muon trigger is 18 GeV, the thresholds for the 
electron trigger are 20 GeV or 22 GeV, depending on 
the data-taking period due to changing LHC luminos¬ 
ity conditions. The event reconstruction makes use of 
kinematic variables such as transverse momentum (px), 
energy in the transverse plane (At) and pseudorapidity 
( 77 ) of photons, leptons (e and p) and jets (j) as well as 6 - 
tagging information, and missing transverse momentum 
(A“ iss ). 

The selected events are required to contain a recon¬ 
structed primary vertex with at least five associated 
tracks, each with px > 0.4 GeV. The primary vertex 
is chosen as the vertex with the highest X)Pt over 
associated tracks. 

Photons are required to have Ax > 20 GeV and p < 
2.37, excluding the transition region between the barrel 
and endcap calorimeters at 1.37< |p| <1.52, and must 
satisfy tight identification criteria [38, 39]. Specifically, 
requirements on the electromagnetic shower shapes [40] 
are applied to suppress the background from hadron de¬ 
cays (e.g. 7T° —► 77 decay leads to two overlapping show¬ 
ers as opposed to a single shower produced by a prompt 
photon). 

Electrons [41] are reconstructed by matching energy 
deposits in the electromagnetic calorimeter with tracks in 
the ID, and are required to have Ax > 25 GeV and |p| < 
2.47, excluding the transition region between the barrel 
and endcap calorimeters. Muons [42] are reconstructed 
by matching tracks in the ID with tracks measured in 
the muon spectrometer, and are required to have px > 
20 GeV and \rj\ < 2.5. 

Leptons are required to be isolated to reduce the num¬ 
ber of lepton candidates that are misidentified hadrons 
or non-prompt leptons. To calculate the isolation of 
electrons in the calorimeter, the Ax deposited in the 


calorimeter in a cone of size A R = 0.2 around the elec¬ 
tron is summed, and the A T due to the electron itself 
is subtracted. The scalar sum of px of tracks with 
px > 1 GeV originating from the primary vertex in 
a cone of A R = 0.3 around the electron direction is 
also measured, excluding the electron track. Selection 
requirements are parameterized as a function of the elec¬ 
tron 77 and Ax and applied to these two isolation variables 
to ensure a constant efficiency of the isolation criteria of 
90% (measured on Z —► e + e“ data) over the entire ( 77 , 
Ax) range. For muons, the transverse energy deposited 
in the calorimeter in a cone of A R = 0.2 around the muon 
direction is required to be less than 4 GeV, after subtrac¬ 
tion of the Ax due to the muon itself. The scalar sum of 
the transverse momenta of tracks in a cone of A R = 0.3 
is required to be less than 2.5 GeV after subtraction of 
the muon track px- The efficiency of the muon isolation 
requirements is of the order of 86 % in simulated tt events 
in the muon+jets channel. 

Jets [43] are reconstructed from topological clus¬ 
ters [44, 45] of energy deposits in the calorimeters us¬ 
ing the anti-fc t [46] algorithm with a distance parameter 
R = 0.4. Jets selected for the analysis are required 
to have px > 25 GeV and |p| < 2.5. In order to re¬ 
duce the background from jets originating from pile-up 
interactions, the jet vertex fraction, defined as the sum 
of the px of tracks associated with the jet and originat¬ 
ing from the primary vertex divided by the sum of the 
px from all tracks associated with the jet, is required 
to be greater than 0.75. Since electrons and photons de¬ 
posit energy in the calorimeter, they can be reconstructed 
as jets. The jet closest to an identified electron in 77 —<^> 
space is rejected if A R(e,j) < 0.2 [47]. Similarly, any 
jet within AR(p/,j) = 0.1 of an identified photon is dis¬ 
carded. To suppress muons from heavy-flavor hadron de¬ 
cays inside jets, muon candidates within A Rip. j) < 0.4 
are rejected [47]. 

Jets containing a 6 -hadron are identified with a 6 - 
tagging algorithm [48-50] using impact parameter and 
vertex position measurements from the inner detector as 
inputs to a neural network; 6 -tagged jets are required to 
satisfy a selection that is 70% efficient for 6 -quark jets in 
simulated tt events. The misidentification rate of light- 
flavor partons (u, d , s-quark or gluon) is in the range 
from 1% to 3%, depending on the jet px and 77 [48]. 

The transverse momentum of the neutrinos produced 
in the top-quark decay chains, measured as missing trans¬ 
verse momentum, is reconstructed from the vector sum of 
the transverse momenta corresponding to all calorimeter 
cell energies contained in topological clusters [43] with 
|r?| < 4.9, projected onto the transverse plane. Contribu¬ 
tions to A“ lss from the calorimeter cells associated with 
physics objects (jets, leptons, photons) are calibrated ac¬ 
cording to the physics object calibration [51]. The con¬ 
tribution to A™ ISS from the px of muons passing the se¬ 
lection requirements is included. Calorimeter cells con¬ 
taining energy deposits above noise and not associated 
with high-px physics objects are also included. 
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Figure 1: Comparison of distributions in data (points) versus expectation (stacked histograms) for the ti selection 
(see text). The electron transverse energy ( E t) in the electron channel is shown on the left, missing transverse 
momentum (I?™ lss ) in the muon channel is shown on the right. The contribution from multijet production and its 
uncertainties are estimated using a data-based technique (see Sec. VIIB). Other contributions are estimated using 
Monte Carlo simulations. The uncertainty band includes statistical and systematic uncertainties. The systematic 
uncertainties include those on lepton, jet, E™ lss , and 6 -tagging modeling, as well as systematic uncertainties on the 
multijet background estimate. The last bin contains any overflow. 
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Figure 2: Distributions for the ti y selection (see text). The photon candidate transverse energy (£r) distribution in 
data (points) is compared to the expectation (stacked histograms) for the electron (left) and muon (right) channels. 
The contribution from multijet +7 production and its uncertainties are estimated using a data-based technique (see 
Sec. VIIB). The remaining contributions are estimated using Monte Carlo simulations. Other backgrounds (labeled 
as ‘Other bck.’) include contributions from Z + jets, dibosons and single-top-quark production. The contribution 
from tt production with prompt photons (labeled as ‘tty’) is estimated using the WHIZARD tty Monte Carlo 
simulation. The contribution from ti events with electrons and hadrons misidentified as prompt photons is obtained 
using inclusive ti Monte Carlo simulation. The uncertainty band includes statistical and systematic uncertainties. 
The systematic uncertainties include those on photon, lepton, jet, E™ lss , and 6 -tagging modeling, as well as 
systematic uncertainties for the multijet background estimate. The last bin contains any overflow. 


Top-quark-pair candidate events are selected by re¬ 
quiring exactly one lepton i (where i is an electron 
or muon) and at least four jets, of which at least 
one must be 6 -tagged. To reduce the background 
from multijet processes, events in the electron chan¬ 
nel are required to have E™ lss > 30 GeV, where E™ lss 
is the magnitude of the missing transverse momen¬ 


tum E™ 1SS , and a IV-boson transverse mass rm(W) > 
35 GeV. Thi g. -jrap'irveisg mass is defined 

as tot(IV) = x E^ iss (l — cos</>), where p ^ is the 

transverse momentum of the lepton and <fi is the az¬ 
imuthal angle between the lepton direction and the miss¬ 
ing transverse momentum vector. Similarly, events in 
the muon channel are required to have E™ lss > 20 GeV 
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and m,T(W) + E™ lss > 60 GeV. Representative distribu¬ 
tions of kinematic variables for this selection are shown 
in Fig. 1. 

The analysis of tty production is performed on the 
subset of selected tt candidate events that contain at 
least one photon candidate. To suppress the contribu¬ 
tions from photons radiated from leptons, photon can¬ 
didates with A< 0.7 are discarded. Events with 
a jet closer than Ah’( 7 , j) = 0.5 in 77 -^ space to any 
photon candidate are discarded, as those photons have a 
reduced identification efficiency. In addition, to suppress 
the contribution from Z{—> e + e _ )+jets production with 
one electron misidentified as a photon, the ey invariant 
mass m ei is required to be | m ei — mz | > 5 GeV, where 
TOz = 91 GeV. This selection yields totals of 140 and 
222 events in data in the electron and muon channels 
respectively. In Fig. 2 the photon candidate E T distri¬ 
butions for this selection are compared to predictions for 
the electron and muon channels. 

Corrections are applied to simulated samples when cal¬ 
culating acceptances to account for observed differences 
between predicted and observed trigger, photon and lep¬ 
ton reconstruction and identification efficiencies and jet 
6 -tagging efficiencies and mistag rates, as well as smear¬ 
ing to match jet [52], photon and lepton energy resolu¬ 
tions in data [42, 53]. 


V. DEFINITION OF THE FIDUCIAL PHASE 
SPACE AND CROSS SECTION 

To allow a comparison of the analysis results to the¬ 
oretical predictions, the cross section measurement is 
made within a fiducial phase space defined in Monte 
Carlo simulation for tty decays in the single-lepton (elec¬ 
tron or muon) final state. The particle-level prediction 
is constructed using final-state particles with a lifetime 
longer than 10 ps. 

Photons are required to originate from a non-hadron 
parent, which is equivalent to the requirement for pho¬ 
tons to originate from a top-quark radiative decay or top- 
quark radiative production. Photons are required to have 
p T > 20 GeV and \r]\ < 2.37. 

Leptons are defined as objects constructed from the 
four-momentum combination of an electron (or muon) 
and all nearby photons in a cone of size A R = 0.1 in rj—ifi 
space centered on the lepton. Leptons are required to 
originate from a non-hadron parent, which is equivalent 
to the requirement for leptons to originate from the t —» 
Wb -A Ivb decays. Leptons are required to have pr > 
20 GeV and \rj\ < 2.5. 

Decays of tty to the dilepton final states, as well as 
decays to the single-lepton final state with an electron or 
muon coming from a r —> tvv T decay are considered as 
non-fiducial and are corrected for when calculating the 
cross section. 

The anti -kt [46] algorithm with a distance parame¬ 
ter R = 0.4 is used to form particle-level jets from 


all particles with a lifetime longer than 10 ps, excluding 
muons and neutrinos. Particles arising from pile-up in¬ 
teractions are not considered. Jets are required to have 
Pt > 25 GeV and \r/\ < 2.5. 

The removal of overlapping particles is performed in 
a manner consistent with the object and event selection 
described in Sec. IV. Any jet with A R(e,j) < 0.2 or 
Ai?(y,j) < 0.1 is discarded; any muon with A R(fi,j) < 
0.4 is discarded. To suppress the contribution of photon 
radiation off a charged lepton, photons within Ai?(y, £) < 
0.7 are discarded. 

For the determination of the tty fiducial cross sec¬ 
tion crJA), exactly one lepton (electron or muon), at least 
one photon, and four or more jets are required. At 
least one jet must match a 6-hadron. All simulated 6- 
hadrons that are generated with pt > 5 GeV are con¬ 
sidered for the matching, and are required to satisfy 
Ai?(6-hadron, j ) < 0.4. Events with A/?(y,j) < 0.5 are 
discarded. 

The fiducial cross section is calculated as 

= N s /(e ■ C). The number of estimated tty signal 
events is N s = N — Nf,, where N and Nb are the number 
of observed tty candidate events in data and the esti¬ 
mated number of background events respectively. The 
efficiency e is determined from tty Monte Carlo simula¬ 
tion as the ratio of the number of all events passing the 
tty event selection to the total number of events gener¬ 
ated in the fiducial region. It is 17.8 ± 0.5 (stat.)% for 
the electron channel and 34.3± 1.0 (stat.)% for the muon 
channel. These numbers include kinematic and geomet¬ 
ric acceptance factors, as well as trigger, reconstruction 
and identification efficiencies. The efficiency values also 
account for migrations into and out of the fiducial phase 
space. 

VI. ANALYSIS STRATEGY 

After the selection more than half of the events do not 
come from tty production. The track-isolation distribu¬ 
tion of the photon candidates is used to discriminate be¬ 
tween signal photons and neutral hadron decays to final 
states with photons and hadrons misidentified as pho¬ 
tons. For simplicity, neutral hadron decays to diphoton 
final states and hadrons misidentified as photons are re¬ 
ferred to hereafter as ‘hadron-fakes’. 

The photon track-isolation variable p'j>° is defined as 
the scalar sum of the transverse momenta of selected 
tracks in a cone of A R = 0.2 around the photon can¬ 
didate. The track selection requires at least six hits in 
the silicon pixel and microstrip detectors, including at 
least one hit in the innermost layer in the pixel detec¬ 
tor (except when the track passes through one of the 
2% of pixel modules known to be not operational), track 
Pt > 1 GeV, longitudinal impact parameter |z 0 | < 1 mm 
and transverse impact parameter |do| < 1 mm computed 
with respect to the primary vertex. The tracks from pho¬ 
ton conversions are excluded. 
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Prompt-photon and background track-isolation tem¬ 
plates are obtained from data as described in Sec. VIB 
and VIC. The total number of events with prompt 
photon-like objects (for simplicity referred to as ‘prompt 
photons’ unless noted otherwise) is extracted using a 
template-based profile likelihood fit. The expected num¬ 
ber of non-tty events with prompt photons, as summa¬ 
rized in Table I, is subtracted to calculate the fiducial 
cross section These steps are incorporated in a like¬ 
lihood fit. 


A. Likelihood description 


bin j spanning the range Vj reads: 


L (plfP | of^, e(6), £, N bl (rii),..., N bn (a n )J = 


N channels 

n 

C=1 


NA 


n 

3=1 


N bkg -syst 

x n 1 ^OLl ) 

1 = 1 


background uncertainties 

X J\f(0 k \0 k ,ae k ) x M(C\C,ac) , 


Poisson expectation 

N syst 


k =1 


luminosity uncertainty 


efficiency/acceptance uncertainties 


(i) 


A binned template fit maximizes the following ex¬ 
tended Poisson likelihood function, representing the Pois¬ 
son probability to observe N data events given an expec¬ 
tation of (N s + N b ) events: 


L (pT | N„ N b ) = (Ns +^ e~ (N ’ +Nb) x P(pT \ N s + N b )x 

n 

n P ( N »i |4)X Peff(e I i) x Plum (£ 1 1). 
i=1 

For a given variable x, P(x\x) is the probability of x 
given x, where x denotes the unconditional maximum- 
likelihood estimate of x. Therefore, P e ff(£|e) describes 
the systematic uncertainties affecting the combined sig¬ 
nal efficiency and acceptance e; Pi um (£ | C) describes the 
uncertainty on the integrated luminosity £; P(N bi \ N bi ) 
describes the uncertainty on the i-tli background com¬ 
ponent bp, n is the number of background sources, N b = 

The modeling of the signal and the different back¬ 
ground sources can be expressed as: 


P{pT\Ns + N b ) = f sb Fs(pT) + (1 - fsb) F *>i (Pt°)i 


where F s (p 1 j,°) and F bi (plj?) are the probability density 
functions (pdf) for the signal and the i- th background 
source respectively, with f sb = N S /(N S + N b ) being the 
signal purity. Each F b . is normalized to the correspond¬ 
ing background expectation N b ./N b . 

Every systematic uncertainty is taken into ac¬ 
count as an independent nuisance parameter mod¬ 
eled by a Gaussian pdf J\f. In the likelihood, 

£ = (^electron channel; ^muon channel) and N bi are Consid¬ 
ered to be functions of the nuisance parameters 9 and d) 
respectively. Taking into account the probability distri¬ 
bution functions modeling the different parameters, the 
expanded form of the likelihood used to fit Wins of the 
distribution for an expectation of Nj events in each 


where Uj is defined as: 


v j = V 3{°ttv N bl (a)), ...,N k (a n )) = 


e c (9)£a?« I dp™F>(p™\<T%) 


fid • 


/ dp^F^\N b ^a t )), (2) 

JVi 


with c ={electron channel, muon channel}, and i = 
1,... ,IV b kg-sy S t and k = 1 ,...,N syst denoting the sys¬ 
tematic uncertainties on the background and the signal 
efficiency/acceptance respectively. The normal pdf, mod¬ 
eling the nuisance parameter x, is denoted by N{x\x, a x ). 
The p'rfP binning is chosen to minimize the statistical un¬ 
certainty. 

Finally, a profile likelihood ratio A s is built [54, 55] by 
considering the cross section as the parameter of interest 
and all other parameters to be nuisance parameters: 


A s (pt 



L(pT\a^,k9),lMd)) 

L(pT\a?^,m,£,N b (a)) 


Here, for a given parameter x , x is the value of x that 
maximizes the likelihood function for a given crf^. The 
numerator thus depends on the conditional likelihood es¬ 
timator of x, and the denominator depends on the max¬ 
imized (unconditional) likelihood estimator. 


B. Prompt-photon template 

The prompt-photon template models the p'rfP distribu¬ 
tion of prompt photons as well as electrons misidentified 
as photons, from ttj and background processes. While 
the same template is used for prompt photons and elec¬ 
trons misidentified as photons, the possible differences 
are covered by alternative templates used to estimate the 
systematic uncertainties as discussed below. 

Since electron and photon track-isolation distributions 
are expected to be very similar, the electron template 
ydat a ,e ex t rac t e d f rom the electron p™° distribution in 
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Z —> e + e candidate data events. The prompt-photon 
template T^ ata is then derived taking into account the 
differences between electron and photon plfP distribu¬ 
tions as well as differences between the Z —> e + e~ and 
ttj event topologies, as photons from tt' y events are less 
isolated than electrons from Z —>■ e + e _ events. To ob¬ 
tain the prompt-photon template, the electron plfP dis¬ 
tribution in Z —> e + e~ candidate data events is cor¬ 
rected using weights ( Wi ) and templates obtained from 
Z —> e + e~ (T s ^’ e ) and tt 7 (T s ^’ 7 ) Monte Carlo simu¬ 
lations in twelve px x V bins (indexed by i ): 


T idata mdata.c , ( mMC.'Y rnMC,6 \ 

sig = T sig + 2^ Wi ) ■ 

i=p r £ ,rf bins 

The three px bins are defined as 20 GeV < px < 
30 GeV, 30 GeV < p T < 50 GeV, p T > 50 GeV. 
The four rj bins are defined as 0.0 < |p| < 0.6, 0.6 < 
|p| < 1.37, 1.52 < \r}\ < 1.81 and 1.81 < |p| < 2.37. 
The relative weight for each bin i is calculated from the 
photon Et and 77 spectra of the ttj Monte Carlo sam¬ 
ple. The prompt-photon template, labeled as ‘Nominal’, 
is shown in Fig. 3. It is shown along with an electron p^° 
template obtained from Z{—> e + e“)+ > 4-jets candidate 
data events, and a prompt-photon plfP template obtained 
directly from tij Monte Carlo simulation. The latter two 
templates are used to estimate systematic uncertainties 
on the measured cross section due to the choice of the 
prompt-photon template. 
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Figure 3: Comparison of the nominal prompt-photon 
track-isolation (px°) template with the template 
obtained from data using a Z(—> e + e _ )+ > 4-jets 
selection, and with the template obtained from tt'y 
simulation. The distributions show the probability 
P(pf°| 7 ) of observing a photon in a given plfP bin per 
GeV. The last bin contains any overflow. 


C. Background template 

Contributions from background sources with non¬ 
prompt photons are described by a single template. 
This background template is extracted from a multijet 
data sample by inverting requirements on photon shower 
shape variables as described in Sec. VIC1. This set of 
events is referred to as the ‘hadron-fake control region’. 
A correction is applied to account for the prompt-photon 
contribution in the background template as described in 
Sec. VIC2. 


1. Derivation 

The hadron-fake control region is obtained from mul¬ 
tijet events that are required to have either at least 
two jets with px > 40 GeV and at least two addi¬ 
tional jets with px > 20 GeV, or at least five jets with 
px > 20 GeV. Non-pronrpt photon candidates are iden¬ 
tified by inverting requirements on the electromagnetic 
shower shapes [40] . The background template shapes are 
determined separately in the four photon 77 bins and three 
photon Et bins defined in Sec. VIB. The photon Et 
distributions are consistent across different 77 regions, so 
?7 and Et: dependencies of the background template are 
treated separately. 

To match the expected px and 77 distributions of non¬ 
prompt photons in the signal region, these seven tem¬ 
plates are weighted using 77 and px distributions of non¬ 
prompt photon candidates in tt candidate events in data. 
The resulting background template (labeled as ‘Nominal 
template T^* a ’) is shown in Fig. 4. 

2. Prompt-photon contribution to the background template 

While the nominal background template is extracted 
using a data-based procedure as described above, the 
prompt-photon contamination in the background tem¬ 
plate is obtained using a combination of data and Monte 
Carlo information. 

Multijet simulation is used to obtain a Monte Carlo 
template modeling the isolation distribution of hadrons 
misidentified as photons, T^ c , by applying the same ob¬ 
ject and event selection as for the nominal background 
template, as described in Sec. VIC1. A subset of the 
events used to construct Tj^ c is selected by the require¬ 
ment that those events do not contain any simulated true 
high-px prompt photons. This subset is used to build 
a template {Tjf G ) which models the isolation distribu¬ 
tion of hadrons misidentified as photons without any true 
prompt-photon contribution. 

Figure 4 shows the comparison of T^ c to the data- 
based background template. The systematic uncertainty 
in each p^° bin of is assigned so that data (Tj^ a ) 
and simulation (Tj^ c ) are in agreement. This uncertainty 
















p®° [GeV] 

Figure 4: A comparison of data-based Tj^ a and 
simulation-based TV C track-isolation background 
templates is shown in the upper panel. The 
distributions show the probability -P(pf°l 7 ) °f observing 
a photon in a given bin per GeV. The ratio of the 
two templates is shown in the lower panel. The hatched 
band shows the total uncertainty. The last bin contains 
any overflow. 

is conservatively taken to be the same for all pifP bins and 
is evaluated to be 27% on values of T^ a (p^°). 

The prompt-photon contamination is then extracted 
from data by maximizing the following extended likeli¬ 
hood function Lf , representing the probability to observe 
N data events in the hadron-fake control region given an 
expectation of n exp : 

n N 

L f = -^re"“ p x 0 [(1 - /) T™ c + /T s t ata ] x 

M{6\0,a e ), (3) 

where T s ^| ta is the prompt-photon template and / is the 

fraction of prompt photons. The parameter 6 is the 
nuisance parameter modeling the systematic uncertainty 
due to the differences between T^‘ a and Tj^ c . The 
fraction of prompt photons is distributed according to 
a Gaussian pdf Af(d\d,<jg) with mean 0=1 and width 
ag = 27%. The result of the fit is shown in Fig. 5, and / 
is determined to be (6.1+J;g) x 10 -2 . The uncertainties 
are obtained at the 68% confidence level (CL) by con¬ 
structing the confidence belt with the Feldman-Cousins 
technique [56] using pseudoexperiments. 

Finally, the signal contamination in the background 
template is included in the general likelihood by means 
of a nuisance parameter af a k e modeling the strength of 
the correction: 

T icorr ( 1 A [Aidata , r .. ,7-. dat.a] 

bkg l 7 I J- bkg C^fake / X 1 s ig 

\ -L C^fake ' J J L -I 



Figure 5: Track-isolation background template 
distribution after maximization of the likelihood Lf 
defined in Eq. 3 (top) and normalized residuals 
(bottom). The markers correspond to the nominal 
hadron background template. The stacked filled 
histograms represent the fraction of prompt photons in 
the hadron-fake control region (obtained as / x X^ 3 ) 
and the fraction of hadron-fakes (obtained from the 
simulation-based template as (1 — /) x Tjj c ) as given 
by the fit. The normalized residuals, shown in the 
bottom plot, are defined as the difference between the 
‘Nominal template’ and the sum of (1 — /) x Tjj c and 
/ x Tf-g* 3 -, divided by the total uncertainty ag. The last 
bin contains any overflow. 


The strength factor Of a ke is constrained to 1 by a Gaus¬ 
sian pdf with width cr a = 28% corresponding to the 
largest of the estimated asymmetric uncertainties on /. 
It is then determined from the general likelihood fit in a 
data-based way. 


VII. PROMPT-PHOTON BACKGROUNDS 

To identify prompt-photon and isolated-electron back¬ 
ground contributions to the events selected in the ttf 
analysis, data-based methods and Monte Carlo simula¬ 
tion are used. These background estimates are summa¬ 
rized in Table I and described below. 


A. Electron misidentifled as a photon 

The contribution from events with an electron misiden- 
tified as a photon is estimated using data by applying the 
e —> 7 misidentification rate to tt + e candidate events. 
The measurement of this misidentification probability 
and cross-checks of the method are described below. 

The sample of events with an electron and a photon 
approximately back-to-back in the transverse plane (in 
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Table I: Estimates of the number of selected events with 
prompt photons, or electrons misidentified as photons, 
from various backgrounds to tt' y production, including 
statistical and systematic uncertainties. 


Background source 

Electron channel 

Muon channel 

e —> q misidentification 

29.4 ± 3.0 

41.5 ± 4.6 

Multijet + q 

1.4 ± 1.2 

1.9 ± 1.1 

W'y + jets 

5.4 ± 1.9 

15.6 ± 4.4 

Single-top-quark + q 

1.8 ± 0.3 

3.8 ± 0.4 

Zq + jets 

2.3 ± 1.6 

4.2 ± 3.1 

Diboson 

0.1 ± 0.1 

0.4 ± 0.1 


c t >) with an electron-photon invariant mass m ei close to 
the -Z-boson mass is dominated by Z — 1 e + e~ decays 
in which one of the electrons radiates a high-Ex photon 
while traversing detector material. The probability for an 
electron to be misidentified as a photon is determined in 
data as a function of the electron transverse momentum 
and pseudorapidity using the eq and e + e~ mass distribu¬ 
tions. One electron (tag) is required to match the single¬ 
electron trigger. Another electromagnetic object (probe), 
an electron or photon, is then required to be present and 
give a di-object mass with the tag close to the E-boson 
mass. The eq and e + e _ mass distributions are fit with 
the sum of a Crystal Ball [57, 58] function (for the signal 
part) and a Gaussian function (for the background part) 
to obtain the numbers of ee and eq pairs, N ee and N ej , 
to which several pairs per event can enter. The proba¬ 
bility of an electron being misidentified as a photon is 
measured in q and px bins as / e _>. 7 = N ei /N ee . 

The nominal selection for the signal ttq region is mod¬ 
ified by replacing the photon requirement by an extra¬ 
electron requirement. This extra electron (ef) must ful¬ 
fill the photon kinematic selection, Ex(ef) > 20 GeV 
and |??( e f)l < 2.37, excluding the transition region 
between the barrel and endcap calorimeters at 1.37< 
|p(ef)| <1.52. To estimate the contribution from an elec¬ 
tron misidentified as a photon, these ‘tt + e’ events are 
reweighted according to the probability of the extra elec¬ 
tron being misidentified as a photon. This procedure 
gives 29.4 ± 3.0 and 41.5 ± 4.6 events in the electron and 
muon channels respectively. 

The misidentification probability is also esti¬ 

mated in Z —> e + e~ Monte Carlo simulation, so that 
a closure test can be performed. The number of back¬ 
ground events in simulation that pass the ttq event 
selection is estimated using generator-level information 
about how the photon is produced. These events are 
weighted with the data-to-simulation correction factors 
Sex = found typically to be within 10 % of 

unity. This estimate is found to be in agreement with 
reweighting the events that pass the ‘tt - 1 - e’ event selec¬ 
tion in Monte Carlo simulation according to f e ~>. y , i.e. 
effectively using the data-based approach in the Monte 
Carlo simulation. 


B. Multijet 4- photon 

The background contribution from multijet events with 
associated prompt-photon production is estimated using 
the data-based matrix method discussed in more detail in 
Ref. [59] . In this method, two sets of lepton selection cri¬ 
teria are defined. The ‘tight’ selection criteria are used to 
identify leptons in ttq candidate events. In the ‘loose’ se¬ 
lection criteria the lepton isolation requirements are dis¬ 
regarded, and looser identification requirements [40] are 
applied for electrons. 

The number of selected ttq candidate events is ex¬ 
pressed as a sum of those with prompt leptons and 
those with ‘fake leptons’ (non-prompt leptons or hadrons 
misidentified as leptons). Identification efficiencies for 
prompt leptons are measured in Z —1 t + t~ (t = e, //) 
data candidate events, whereas the efficiency for fake 
leptons to be identified as ‘tight’ leptons is measured in 
a multijet data sample. The number of ttq candidate 
events with at least one non-prompt lepton candidate is 
estimated using this information [59]. 

A template fit to the photon plfP distribution is used to 
determine the prompt-photon fraction in selected ‘mul¬ 
tijet + q’ events. The ‘multijet + q’ event selection 
is similar to the ttq selection except that ‘loose’ lepton 
identification criteria are used instead of the ‘tight’ cri¬ 
teria. Assuming that the prompt-photon fraction does 
not depend on the lepton identification criteria (‘loose’ or 
‘tight’), this prompt-photon fraction is then used to esti¬ 
mate the contribution of the multijet + prompt-photon 
process to the ttq event selection. This results in 1.4 ± 
1.2 and 1.9 ± 1.1 events expected for the electron and 
muon channels respectively. 


C. Wq + jets production 

Background from kEq+jets production is estimated by 
extrapolating the number of IVq+jets candidate events in 
data from a control region (CR) to the ttq signal region 
(SR) using fVq + jets Monte Carlo simulation [60]. In 
the control region the lepton, photon, E™ lss and mx(W) 
selection criteria are the same as in the nominal ttq selec¬ 
tion. To enrich the control region in Wq+jets, events are 
required to have one, two or three jets, and a 6 -tagging 
veto is applied. 

To estimate the prompt-photon contribution, it is as¬ 
sumed that the fraction of prompt photons is the same 
in the CR and SR. To verify this assumption, a template 
fit to the photon distribution is performed, and the 
prompt-photon fraction in data and simulation is found 
to be independent of the jet multiplicity. 

To suppress the Z + jets background contribution 
in the CR, the m e7 requirement is extended to 
|m e7 — rriz \ > 15 GeV. The multijet + q contribution 
to the Wq+jets background in the CR is estimated us¬ 
ing the matrix method as described in Sec. VIIB. The 
number of IVq+jets events with prompt photons in the 
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CR is estimated using a template fit to the photon 
distribution. 

Other contributions to the LVy+jets CR are estimated 
using simulation, where events are separated into two 
classes, one with a prompt photon, the other with an 
electron misidentified as a photon. To obtain the e —> y 
contribution, the s e7 correction factors (Sec. VII A) are 
used. A comparison of data and expectation in the CR 
is presented in Table II. 

Table II: Data and simulated background yields in the 
LVy+jets data control region. The number of events 
with a prompt photon in data (labeled as ‘Events with 
prompt 7 ’) is estimated from the total number of 
lV 7 + jets candidate events in the control region 
(labeled as ‘IVy + jets control region’) using template 
fits. Background yields are estimated using Monte Carlo 
(MC) simulation, except for the multijet + 7 yield. The 
resulting number of IT 7 candidate data events, as well 
as the MC prediction for the number of IVy events are 
shown. To obtain the LVy + jets background to the tty 
selection, the number of IVy candidate data events is 
extrapolated into the signal region using Monte Carlo 
simulation. The uncertainties include both the 
statistical and systematic uncertainties. 



Electron channel 

Muon channel 

LVy+jets control region 

3410 

8394 

Events with prompt 7 

2412 

5540 

tty 

82 ± 16 

161 ± 32 

Z + jets 

160 ± 90 

620 ± 330 

Diboson 

13 ± 3 

26 ± 7 

Single-top-quark 

9 ± 2 

20 ± 5 

e —► 7 misidentification 

380 ± 110 

330 ± 40 

Multijet + 7 

60 ± 30 

350 ± 70 

Total background 

700 ± 140 

1510 ± 340 

LVy estimate 

1710 ± 180 

4030 ± 390 

IT 7 MC expectation 

1860 ± 200 

3930 ± 390 


The number of LVy+jets candidate events in the CR 
(< 3 jets) is extrapolated to the jet multiplicity of the 
SR, > 4 jets [59]. To extrapolate from the LVy+jets 
event selection, which has a 5-tagging veto, to the SR, the 
heavy-flavor quark content is studied in data in events 
with a W boson and two jets. The heavy-flavor quark 
content is then extrapolated from the IVy + 2-jets re¬ 
gion into the SR using the LVy + jets simulation [59, 60]. 
This extrapolation accounts for the difference in flavor 
composition between the IVy + 2-jet and LVy+ > 4-jet 
samples as well as for differences in the per-flavor event 
tagging probabilities, which may lead to different event 
rates after 6 -tagging. The LVy+jets background estimate 
is 5.4 ± 1.9 and 15.6 ± 4.4 events for the electron and 
muon channels respectively. 

Monte Carlo modeling uncertainties in the estimate of 
the background from IVy + jets production include con¬ 
tributions from the estimated number of events with elec¬ 
trons misidentified as photons (which is known to 10 %) 


and from cross section uncertainties (e.g. a 48% uncer¬ 
tainty for Z +jets contributions, which corresponds to the 
error on the normalization of Z+ jets in the four-jet bin 
from the Berends-Giele scaling [60]). 

D. Other background sources 

The single-top-quark, Z+ jets, and diboson contribu¬ 
tions are estimated from simulation and normalized to 
theoretical calculations of the inclusive cross sections. 

The single-top-quark production cross section is nor¬ 
malized to the NLO+NNLL prediction: the i-channel to 
64.6^7 pb [61], the s-channel to 4.6±0.2pb [62], and the 
TVi-channel to 15.7±1.2pb [63]. The Z+ jets background 
is normalized to the NNLO QCD calculation for inclusive 
Z production [64] and the diboson background is normal¬ 
ized to the NLO QCD cross section prediction [65]. 


VIII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties may affect the shapes of the 
prompt-photon and background templates, the esti¬ 
mates of background components with prompt photons 
and with electrons misidentified as photons, as well as 
the efficiencies, acceptance factors and the luminosity. 

The total effect of each systematic uncertainty on the 
cross section is evaluated using ensemble tests. For each 
systematic uncertainty i. pseudodata are generated from 
the full likelihood while keeping all parameters fixed to 
their nominal values except for the nuisance parameter 
corresponding to the systematic uncertainty source. For 
each set of pseudodata, a template fit is performed allow¬ 
ing all parameters of the likelihood (nuisance parameters, 
signal cross section) to vary. The distribution of cross 
sections obtained form a Gaussian pdf with a width that 
gives the uncertainty in the cross section due to the i -th 
systematic uncertainty. This method provides an esti¬ 
mate of the effect of each uncertainty on the cross section 
as shown in Table III. Uncertainties obtained with this 
method are by construction symmetric. All systematic 
uncertainties are described in the following. 

Table III: Summary of systematic uncertainties on the 
tty fiducial cross section, 


Uncertainty source 

Uncertainty [%] 

Background template shapes 

3.7 

Signal template shapes 

6.6 

Signal modeling 

8.4 

Photon modeling 

8.8 

Lepton modeling 

2.5 

Jet modeling 

16.6 

6 -tagging 

8.2 

modeling 

0.9 

Luminosity 

1.8 

Background contributions 

7.7 
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A. Template shapes 

The contribution to the systematic uncertainty on crM 
due to the template shape modeling amounts to 7.6% in 
total. Of this, the background template shape modeling 
uncertainty amounts to 3.7% of the cross section, and the 
prompt-photon template uncertainty amounts to 6.6%. 

The prompt-photon template shape systematic uncer¬ 
tainty is estimated with pseudoexperiments by replac¬ 
ing the nominal prompt-photon template with alterna¬ 
tive templates shown in Fig. 3: (a) an electron pifP tem¬ 
plate obtained from Z(—> e + e“)+ > 4-jets candidate 
data events (4.1% systematic uncertainty is obtained) 
and (b) a prompt-photon p™ 0 template obtained directly 
from tty Monte Carlo simulation (6.6% systematic uncer¬ 
tainty is obtained). The larger of the two uncertainties 
is used as the systematic uncertainty. 

The systematic uncertainty associated with the 
reweighting of the background template is estimated by 
varying within their uncertainties the non-prompt pho¬ 
ton px- and ^-distributions that are used for reweighting. 
The effect of this systematic uncertainty on the cross sec¬ 
tion measurement is found to be negligible. To estimate 
the systematic uncertainty due to the amount of prompt- 
photon contamination in the background template (as de¬ 
scribed in Sec. VIC), the corresponding nuisance param¬ 
eter afake is sampled using a Gaussian pdf with a width of 
<r arake = 28% corresponding to its estimated uncertainty. 
The systematic uncertainty on the cross section is esti¬ 
mated to be 3.7%. All template-shapes uncertainties are 
taken as fully correlated between the electron channel 
and the muon channel. 


B. Signal modeling 

The uncertainty on the tty cross section (as defined in 
Sec. V) due to the modeling of the signal is estimated 
to be 8.4%. The estimate is obtained by varying the se¬ 
lection efficiency with respect to the nominal tty Monte 
Carlo sample which includes event migrations into and 
out of the fiducial region. This uncertainty includes a 
comparison of MadGraph with WHIZARD (1.7%), as well as 
a comparison of the MadGraph tty samples with different 
QED FSR settings (3.4%) as explained in Sec. III. The 
renormalization and factorization scales are also varied, 
leading to an uncertainty of 1.1%. To assess the effect 
of different parton shower models, predictions from the 
MadGraph+HERWIG sample are compared to predictions 
from the MadGraph+PYTHIA sample, leading to an uncer¬ 
tainty of 7.3%. In addition, studies of ti samples with 
varied color reconnection (0.2%) and underlying event 
(0.9%) settings lead to small contributions. The un¬ 
certainty associated with the choice of the CTEQ6L1 
PDF set is evaluated from an envelope of calculations 
using the PDF4LHC prescription [66] by reweighting the 
CTEQ6L1 LO PDF used in the generation of the tty 
WHIZARD sample with MSTW2008 [67, 68], CT10 [34, 69] 


and NNPDF2.0 [70] NLO PDF sets and amounts to 1.1%. 
All signal-modeling uncertainties are taken as fully cor¬ 
related between the electron channel and the muon chan¬ 
nel. 


C. Detector modeling 

The systematic uncertainty on the cross section due to 
photon modeling is 8.8%. It is estimated from the pho¬ 
ton identification (7.3%) [38], the electromagnetic energy 
scale (2.7%) and the resolution (4.0%) systematic uncer¬ 
tainties [53]. 

The systematic uncertainty on the cross section due to 
lepton modeling is 2.5%. It is estimated separately for 
the electron and muon channels from the lepton trigger 
(0.3% and 1.7%), reconstruction (0.5% and 0.4%) and 
identification (1.2% and 1.0%) efficiency uncertainties, as 
well as from those on the energy scale (0.3% and 0.3%) 
and resolution (0.1% and 0.7%). 

The systematic uncertainty on the cross section due 
to jet modeling is 16.6%. It is estimated taking into 
account the following contributions. The largest ef¬ 
fect comes from the energy scale (15.0%) uncertainty 
which is estimated by combining information from the 
single-hadron response measured with in-situ techniques 
and with single-pion test-beam measurements [52]. The 
jet energy resolution (6.5%) uncertainty is estimated 
by smearing the jets in simulation by the uncertainty 
as measured with the dijet balance and bisector tech¬ 
niques [71]. The uncertainty on jet reconstruction effi¬ 
ciency (1.0%), which is defined relative to jets built from 
tracks reconstructed with the ID, is also considered [43]. 
The jet vertex fraction uncertainty is found to be 2.6%. 

The systematic uncertainty on the cross section due 
to 6-tagging modeling is 8.2%. It is dominated by the 
contribution due to the efficiency (8.1%) [49] with a small 
contribution due to the mistag probability (1.1%) [48]. 

Systematic uncertainties on the energy scale and res¬ 
olution of leptons, jets and photons are propagated to 
Jf“ lss . Additional I?™ lss uncertainties [51] also taken into 
account are contributions from low-px jets and from en¬ 
ergy in calorimeter cells that are not included in the re¬ 
constructed objects (0.3%), as well as any dependence on 
pile-up (0.9%). 

All detector-modeling systematic uncertainties except 
for the lepton-modeling uncertainties are taken as fully 
correlated between the electron channel and the muon 
channel. The lepton-modeling uncertainties are taken as 
uncorrelated between channels. 

The effect of the luminosity uncertainty on the cross 
section amounts to 1.8% [4]. 

D. Background contributions 

The total systematic uncertainty originating from the 
non-tty background contributions with prompt photons 
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Figure 6 : Results of the combined likelihood fit using the track-isolation ( pif?) distributions as the discriminating 
variable for the electron (left) and muon (right) channels. The contribution from tt'y events is labeled as ‘Signal’, 
prompt-photon background is labeled ‘7 background’, the contribution from hadrons misidentified as photons (as 

estimated by the template fit) is labeled as ‘Hadron fakes’. 


or electrons misidentified as photons is estimated to be 
7.7%. This uncertainty includes the following: electrons 
misidentified as photons (5.0%), kFq+jets (5.4%), as well 
as multijet + photon (1.5%), Z 7 +jets (1.3%), diboson 
(0.4%) and single-top-quark (0.4%) processes. The vari¬ 
ous sources of uncertainty on the background estimates 
quoted above are described in the following paragraphs. 

For background estimates obtained using simulation, 
uncertainties on the cross section predictions are taken 
into account. Cross section systematic uncertainties are 
considered as fully correlated between the electron and 
the muon channels. However, the corresponding statis¬ 
tical uncertainty is taken as uncorrelated. For Z 7 + jets, 
single-top-quark and diboson contributions the cross sec¬ 
tion systematic uncertainty is negligible with respect to 
the statistical uncertainty. 

The systematic uncertainty on the probability of an 
electron to be misidentified as a photon as described in 
Sec. VIIA is obtained by varying the fit functions and 
the ee and eq mass windows in Z —> e + e _ candidate 
events in data. This uncertainty is estimated to be about 
10 % of the background estimate and it is taken as fully 
correlated between the electron channel and the muon 
channel. 

For the multijet + photon background described in 
Sec. VIIB, the uncertainty is about 90% for the electron 
channel and 60% for the muon channel. It is dominated 
by the statistical uncertainty due to the small number of 
events in the data samples and the systematic uncertain¬ 
ties on the matrix method (50% for the electron channel 
and 20% for the muon channel) [59] . Those uncertainties 
are taken as uncorrelated between the two channels. 

The systematic uncertainties on the kFq+jets back¬ 
ground are dominated by the extrapolation from the con¬ 
trol region (dominated by kFq+jets) to the signal region 
due to different event topologies in the two regions in 


terms of the total number of jets and the number of 
heavy-flavor jets. The uncertainties due to the extrap¬ 
olation are 27% in the electron channel and 23% in the 
muon channel and are dominated by the uncertainty on 
the knowledge of the flavor compositions of the kF+jets 
events and the overall kF+jets normalization for differ¬ 
ent jet multiplicities [59, 60]. Those uncertainties are 
taken as fully correlated between the electron channel 
and the muon channel. The statistical uncertainty on 
the number of events in the kFq+jets control region is 
taken as uncorrelated between the two channels. System¬ 
atic uncertainties on the multijet+photon contribution to 
the kFq+jets event selection, as well as uncertainties on 
Monte Carlo modeling of tt, Z+ jets, diboson, and single- 
top-quark processes are taken into account [47]. 


IX. RESULTS 

Totals of 140 and 222 ttq candidate data events are 
observed in the electron and muon channels respectively. 
The numbers of background events extracted from the 
combined likelihood fit are 79 ± 26 for the electron 
channel and 120 ± 39 for the muon channel. The num¬ 
bers of tij signal events are determined to be 52 ± 14 
and 100 ± 28. The results include statistical and sys¬ 
tematic uncertainties. These numbers are summarized in 
Table IV, and the pif° distributions are shown in Fig. 6 . 

Using the asymptotic properties [72] of the likelihood 
model, the test statistic for the no-signal hypothesis is ex¬ 
trapolated to the likelihood ratio value observed in data 
(14.1) to determine the p-value of pg bs = 5.73 x 10 -8 . 
The process tty in the lepton-plus-jets final state is ob¬ 
served with a significance of 5 . 3 <r away from the no-signal 
hypothesis. 

The tt') fiducial cross section together with its total 
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Table IV: Number of tt -y signal and background events 
extracted from the likelihood fit, which is performed for 
the electron and muon channels simultaneously. The 
uncertainties are statistical and systematic. The total 
number of tty candidate events observed in data is also 
shown. 


Contribution 

Electron chan. 

Muon chan. 

Total 

Signal 

52 ± 14 

100 ± 28 

152 ± 31 

Hadrons 

38 ± 26 

55 ± 38 

93 ± 46 

Prompt photons 

41 ± 5 

65 ± 9 

106 ± 10 

Total background 

79 ± 26 

120 ± 39 

199 ± 47 

Total 

131 ± 30 

220 ± 48 

351 ± 59 

Data candidates 

140 

222 

362 


uncertainty is obtained from the profile likelihood ra¬ 
tio fit to be 63+)® fb. The total systematic uncertainty 
is extracted from y/ (a- syst ® s tat ) 2 - af tat - a\ fb, 

where ac is the luminosity uncertainty; <r s t a t is the pure 
statistical uncertainty, evaluated from the profile likeli¬ 
hood without including nuisance parameters; rx S yst©stat 
is the total uncertainty extracted from the 68 % CL of 
the profile likelihood fit (including nuisance parameters), 
as shown in Fig. 7. 


single-muon final state. Good agreement is found with 
the predicted cross sections [18, 73] of 48 ± 10 fb and 
47 ± 10 fb obtained from the WHIZARD and MadGraph 
Monte Carlo generators respectively and then normal¬ 
ized by the corresponding NLO/LO AT-factors. In ad¬ 
dition, the cross section measurements are performed 
separately in the electron and muon channels and give 

x BR = 76+]J(stat.) (syst.) ± l(lumi.)fb and 

x BR = 55 ^ 9 °(stat.) In (syst-) ± l(lumi.) fb respec¬ 
tively. 


X. SUMMARY 

The production of tty final states with a photon 
with transverse energy greater than 20 GeV is observed 
with a significance of 5.3 ct in proton-proton collisions at 
yfs = 7 TeV using the ATLAS detector at the CERN 
LHC. The dataset used corresponds to an integrated lu¬ 
minosity of 4.59 fb -1 . The tty cross section per lepton 
flavor, determined in a fiducial kinematic region within 
the ATLAS acceptance defined in Sec. V, is measured to 
be af^ x BR = 63 ± 8 (stat.) +*3 (syst.) ± 1 (hum.) fb in 
good agreement with the theoretical prediction. 



Figure 7: Negative logarithm of the profile likelihood as 
a function of the tty fiducial cross section crM x BR 
with (solid line) and without (dashed line) free nuisance 
parameters associated with the systematic 
uncertainties. The horizontal dotted line corresponds to 


a value of — log 




tt^f J 


= 0.5. Intersections of 


this line with the solid (dashed) curve give the ± 1 <t 
total (statistical only) uncertainty interval to the 
measured fiducial tty cross section. 


The tty fiducial cross section times BR per lep¬ 
ton flavor, as defined in Sec. V, is determined to be 
x BR = 63 ± 8 (stat.) (syst.) ± 1 (lumi.) fb, where 
BR is the tty branching ratio in the single-electron or 
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Appendices 

Appendix A: tty Monte Carlo samples 

Signal tt 7 events with single-lepton {lu^qq'bby, 
£ = e, g, t) or dilepton {£v^£'upbby, i/i' = e, g, r) final 
states are simulated with two independent leading-order 
(LO) matrix element (ME) Monte Carlo generators, 
WHIZARD vl.93 [9, 10] and MadGraph v5.1.5.12 [11], 
both using the CTEQ6L1 [12] LO parton distribution 
function (PDF) set. Both calculations take into ac¬ 
count interference effects between radiative top-quark 
production and decay processes. 

1. Leading-order calculations: WHIZARD and MadGraph 

In the WHIZARD tty sample, the minimum transverse 
momentum of all outgoing partons except for the photon 
is set to 10 GeV. The transverse momentum of the pho¬ 
ton is required to be larger than 8 GeV. The invariant 
mass of the photon and any charged particle (u-, d-, c- 
and s-quarks, electrons, muons, and r leptons) is required 
to be larger than 5 GeV. To avoid infrared and collinear 
divergences, the following invariant masses are also re¬ 
quired to be larger than 5 GeV: m(qi,q 2 ), m(gi,qi ), 
m(g 1 ,q 2 ), m(g 2 ,qi), and m{g 2 ,q 2 ), where qi and q 2 are 
the quarks from the hadronic decay of one W boson, and 
g i and g 2 are the gluons initiating the gg —t tty process. 
For each incoming quark Qi (u-, d-, c-, s- and 5-quark), 
the invariant mass to(Q,;, qj ) is required to be larger than 
5 GeV if qj is the same type of parton as Qi. The renor¬ 
malization scale is set to 2 m t , and the factorization scale 
is set to the partonic center-of-mass energy y/§. The cross 
section is 648 fb when summing over all three lepton fla¬ 
vors for the single-lepton (e, g, r) and 188 fb for the 
dilepton tty final states. 

In the MadGraph tty sample, the minimum transverse 
momentum is set to 15 GeV for it-, d-, c- and s-quarks, 
as well as for photons, electrons, muons and r leptons. 
The distance in rj -(f) space between all these particles is 
required to be A R > 0.2. For 5-quarks, no requirement 
is placed on the transverse momentum or on the pseu¬ 
dorapidity. Leptons and photons are required to have 
\r]\ < 2 . 8 , while it-, d-, c- and s-quarks are required to 
have |ry| < 5.0. The renormalization and factorization 
scales are set to m*. The cross section is 445 fb when 
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summing over all three lepton flavors for the single-lepton 
and 131 fb for the dilepton tty final states. 


2. Next-to-leading-order calculation 

The NLO QCD calculation of top-quark pair produc¬ 
tion in association with a hard photon is described in 
Ref. [73] for ^/s = 14 TeV. A dedicated calculation at 
y/s = 7 TeV both at LO and at NLO has been per¬ 
formed for this analysis [18] for the pip -A bg + u ll bjjy 
channel using the same settings for the renormalization 
and factorization scale as in the WHIZARD tty calculation. 

The following NLO input parameters are used: top- 
quark mass mt = 172 GeV, top-quark width T t = 
1.3237 GeV, W-boson mass m w = 80.419 GeV, W- 
boson width = 2.14 GeV, fine-structure constant 
a = 1/137. The strong-coupling constant a s (g) is evalu¬ 
ated using the two-loop running from a s (mz ) as specified 
in the MSTW2008 NLO PDF. Jets are defined using the 
anti-A:* algorithm with a distance parameter R = 0.4. 
The photon is required to be separated from hadronic 
activity as defined in Ref. [74], 

The phase-space requirements used in the yfs = 7 TeV 
theory LO and NLO calculations are described below. 
The muon is required to have p r r(g) > 20 GeV and 
\g{g)\ < 2.5. The missing transverse momentum is re¬ 
quired to be F™ iss > 25 GeV and F“ iss +m^ >60 GeV, 
where to is the IV-boson transverse mass. Jets are 
required to have pr(j) > 25 GeV and \r](j)\ < 2.5. 
The photon is required to have pt(>) > 15 GeV and 
\r](y)\ < 1.37 or 1.52 < \rj{y)\ < 2.37. The objects 
are required to be separated in A R: Ai?(jets) > 0.4, 
Ai?(/x,jets) > 0.4, AR(y,g) > 0.4, Ai?( 7 ,jets) > 0.5. 
The event is required to have -/Vj ets > 4. 

With the above setup and assuming 100% efficiencies, 
ct nl° = 24 . 5 + 8-6 pb and ct lo = 14.7+5-8 pk Uppe r and 

lower values correspond to scale variations by a factor of 
two around g = m t . Therefore, for g = m t the NLO/LO 
AT-factor is 1.67. Similarly, for the WHIZARD Monte Carlo 
sample scales and NLO calculation at the scale of g = mt, 
the NLO/LO I\ -factor is 2.53. 

The LO cross sections calculated with the WHIZARD and 
MadGraph Monte Carlo generators are multiplied by the 
corresponding A-factors in order to compare with the 
experimental measurement. 
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